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Abstract 

We propose a new construction of two-dimensional natural bi-Hamiltonian systems 
associated with a very simple Lie algebra. The presented construction allows us to 
distinguish three families of super-integrable monomial potentials for which one addi- 
tional first integral is quadratic, and the second one can be of arbitrarily high degree 
with respect to the momenta. Many integrable systems with additional integrals of 
degree greater than two in momenta are given. Moreover, an example of a super- 
integrable system with first integrals of degree two, four and six in the momenta is 
found. 



1 Introduction 

The main aim of this paper is to study natural integrable systems 

H^ = 2piP2 + V{q^,q2), (1.1) 
with two degrees of freedom associated with the various representations of Lie algebra 

[N,fl+] = Ki fl+, [N,fl_] = — ?C2fl-, [fl+,fl_] =0 (1.2) 
labelled by parameters Ki,K2 G M, and having the Casimir element 

An obvious representation of this algebra is the algebra of smooth function defined on the 
phase space generated by 

fl+ = qi, fl_ = q2, N = -Kiqipi + K2(j2P2, (1-3) 

with bracket [a,b] := {a,b}, where {■, •} denotes the canonical Poisson bracket in M^. 
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We met with this algebra investigating bi-Hamiltonian structures of two-dimensional 
natural Hamiltonian systems with homogeneous potentials and Newton's equations with 
homogeneous velocity independent forces. Namely, let us consider a natural system given 
by the Hamilton function of the form (|LT) with a monomial potential. More precisely, let 
us consider system governed by Hamiltonian of the form 

H = a+=2prp2 + q^qi (1.4) 

where a and b are real parameters. For such a system the function 

N = apiqi + ^p2qi, (1-5) 

satisfies the following equality 

{N,fl+} = (a + j6)fl+, 
provided that a. and f> are fulfils equation 

fla + bj6 + a + j6 = 0. 

If we can find the remaining generator fl_ = H2, then we have integrable in the Liouville 
sense dynamical system with homogeneous potential. In order to get this additional first 
integral we use the machinery of bi-Hamiltonian geometry, see II221I26II . 

It is worth to notice that the existence of function N with the prescribed property is 
related with a certain symmetry of the system. Let Xf denote the Hamiltonian vector 
field generated by f . Assume that for Hamiltonian \\A) there exists function N such that 
{N, H} = cH, where c is a constant. Then the Hamiltonian vector field is a master 
symmetry of Xh, as [X]v,Xh] = cXu- Sometimes a master symmetry is called a conformal 
symmetry, see [3l. 

Let we assume that X^ with N of the following form 

N := A{qi,q2)pi + B{qi,q2)p2, (1-6) 

where A{q^,q2) and B{qi,q2) are differentiable functions, is a conformal symmetry of X^. 
Then it is easy to show that necessarily A = ccqi and B = ^q2. Thus, equality {N, H} = cH 
implies that 

c = a + j6, afl + + c = 0. (1.7) 

So, we recovered the assumed form of N, see (|1.5|) . Now, the last relation in \1.2\ tells that 
we need an additional integral H2 = fl-. Moreover, we require that {N, H2} = dH2, for a 
certain d G R, and this is equivalent that Xn is also a conformal symmetry of Xhj. 

We are going to work with systems (11.411 with two degrees of freedom because such sys- 
tems appeared as subsystems on invariant manifold of n-dimensional Hamiltonian systems 

The plan of this paper is following. In the next section we start with short descrip- 
tion how Hamiltonian systems of the form (ll.4|) appeared in our investigations of the 
integrability of natural Hamiltonian systems with homogeneous potentials. In Section |3] 
bi-Hamiltonian irregular Poisson manifolds as well as their application for construction of 
first integrals of considered systems are presented. The remaining sections contain results 
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of the integr ability analysis. In Section |4] four families of integrable systems with additional 
first integrals quadratic in the momenta are given. In the next four sections super-integrable 
cases in these families are distinguished. In Section |9] examples of integrable systems with 
additional first integrals of degree greater than two in the momenta are given. In Appendix 
we collected basic facts concerning the Gauss hypergeometric equation which are used in 
this paper. 

2 Monomial potentials 

Let us consider Hamiltonian systems with n degrees of freedom with the Hamiltonian of 
the following classical form 

H = ^EpKn^), (2.1) 

where q = {qi, . . . , and p = {pi, . . . , p„) are the canonical coordinates and momenta; 
the potential V{q) is a homogeneous function of degree /c G Q. The strongest necessary 
integrability conditions for such systems with k ^ Z, were obtained thanks to an application 
of differential Galois methods, see lll5l[T6l[TT1l . To derive conditions of this type one needs 
to know a particular non-equilibrium solution of the considered systems. For the systems 
given by Hamiltonian (12.111 a particular solution can be find in systematic way. Namely, if 
a non-zero d G C" satisfies 

V'{d) = 7d, (2.2) 
for some 7 G C, then the system has a particular solution of the following form 

q{t) = f{t)d, p{t) = q>{t)d, (2.3) 

where q){t) is a scalar function which is defined in the following way. If vector d satisfies 
equation (I2.2D with 7 = 0, then (p(t) := t, and d is called an improper Darboux point of 
potential V. In this case, if G N, and the system is integrable in the Liouville sense, then 
all the eigenvalues of the Hessian matrix V{d) vanish, see lfT9ll . 

If the considered d satisfies equation (|2.2|) with 7 7^ 0, then (p{t) is a solution of equation 
cp = —q)^^^. In this case d is called a proper Darboux point of V, and if the system is 
integrable, then for a given k ^ Z \ {—2,2}, all eigenvalues of V"{d) belong to a certain 
infinite subset of rational numbers Q, see IfTSllTSllTSll . 

As we can see, in spite of the fact that the differential Galois theory is quite involved 
the final result has the form of simple arithmetic restrictions on the eigenvalues of ma- 
trix V"{d). Moreover, for polynomial potentials some relations between eigenvalues of the 
Hessian calculated at all proper Darboux points exist. These relations together with arith- 
metic restrictions on the eigenvalues forced by necessary integrability conditions enable to 
find effectively explicit forms of integrable potentials at least for small n and k. Such a 
systematic analysis was initiated for n = 2 in IfTOl [TTll and later on it was developed for 
n > 2, see lfT9lljU]l . 

However, for some classes of potentials the above approach does not work. In particular 
it is a case if the consider potential does not have any proper Darboux point, and moreover, 
the Hessian matrix V"{d) at each improper Darboux point is nilpotent. For n = 2 an almost 
complete characterisation of such polynomial potentials is given by the following lemma. 
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Lemma 2.1. If a polynomial potential V of degree k > 2 does not have any proper Darboux point, 
then it is either equivalent to the following one 

Vk,i = ci{q2-iqif-\q2 + iqiy, aGC*:=C\{0}. (2.4) 

for some I = 2, . . . ,k — 2, or k = 2s and V has factor {qx ± iq^) with multiplicity s. 

We say that potential V is equivalent to W iff V{q) = W{Aq), where A is n x n matrix 
satisfying AA^ = j6Id„ for a certain jf3 G C*. 

Let us notice that potential V/^j with 2 < I < k — 2 has two improper Darboux points 
di = (l,i) and ^2 = (1, — i). Moreover, V"{di) and V"{d2) are nilpotent. In other words, 
for these potentials we do not have any obstacles for the integrability. 

The potentials of the form (|2.4p look more attractive if we introduce new canonical 
variables 

1 1 

z-i = qi+iq2, Z2 = qi-iq2, yi = ^{pi - ipi), yi = :^{pi + ipi) ■ (2.5) 

In these variables the Hamiltonian function has the form 

H = 2yiy2 + z[zl-'. (2.6) 

Clearly, the systems governed by Hamiltonian function of the form (I1.4D are just generali- 
sation of (|2.6|l for cases when k and / are rational. 



3 Irregular Poisson manifolds 

A bi-Hamiltonian manifold M is a smooth manifold endowed with a pair of compatible 
Poisson bi-vectors P and P' such that 

[[P,P1=0, IP',P1=0, (3.1) 

where ||-, -I is the Schouten bracket. If dimM = 2n, and P is invertible Poisson bivector 
on M, the Nijenhuis operator which is called also the hereditary, or recursion operator, is 
defined as 

Af = P'p-\ 

We say that a bi-Hamiltonian manifold M is regular iff dimM = 2n, P is invertible 
Poisson bivector, and the recursion operator A/" has, at every point, n distinct functionally 
independent eigenvalues. If the recursion operator N" does not have this property then we 
say that bi-Hamiltonian manifold M is irregidar. 

For a regular bi-Hamiltonian manifold M functions 

Hfc = -'-TrA/''^ for k e N, (3.2) 

form a bi-Hamiltonian hierarchy on M, i.e., the Lenard relations hold 

P'dH^ = PdHfc+i , for all k>l. (3.3) 

Generally it is unknown if it is possible to construct in a systematic way a full involutive 
set of functions on irregular bi-Hamiltonian manifolds. 

Anyway, we are going to apply the algebra (|1.2|) to construct integrable systems on 
two-dimensional irregular bi-Hamiltonian manifolds. 
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3.1 The rational Calogero-Moser system 

In order to justify better some aspects of our construction we consider the n-particle rational 
Calogero-Moser model associated with the root system An- It is defined by the Hamilton 
function 



H 



! = 1 



(3.4) 



where a is a coupling constant. The canonical variables {(j,p) satisfy the standard Poisson 
bracket relations {qi, pj} = Sij, and the associated Poisson bivector will be denoted by P. 
This system admits independent and conunuting first integrals 

H^ = /-TrL^ for k = l,...,n, 



where L is the standard Lax matrix 



Pi 



P2 



Vn 



Besides these n integrals of motion the rational Calogero-Moser system admits (n — 1) 
additional functionally independent integrals of motion K„i given by 



= mgiH,n - gmtii, m = 2, n, 



where 



Si 



Y^ql,Hj\, for j = l,...,n. 



! = 1 



According to |[26ll integrals of motion Hj^ and Km can be obtained from the Hamilton 
function H = H2 given by (13.411 as polynomial in momenta solutions of the following 
equations 

PdH = k-^P'd\nHk = {m - 1)-^ P' dlnK^, (3.5) 
where P' is the Poisson bivector compatible with P given by 



P' 



R n 

-H 



+ 



A' 

-A M 



and n X n matrices Yl, A, R and M have the following entries 

n r.1 



H, 



A,- 



and 



R 



dpj 



'1^ (dAu dAkj 



^1 V 



Pk- 
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In this case recursion operator J\f — P'P ^ generates only the Hamilton function 

TrAf'' = 2{2Hf such that PdH = P'dlnH, (3.6) 

and instead of the standard Lenard relations (|3.3|l we have equations l|3.5|). Thus, the de- 
scribed Calogero-Moser system is a super-integrable system on an irregular bi-Hamiltonian 
manifold R^'\ 



3.2 Algebraic construction of two-dimensional bi-Hamiltonian systems on ir- 
regular manifolds 

In order to get integrable systems associated with the Lie algebra (|1.2|) we start with the 
canonical Poisson bivector 









1 


0" 











1 


-1 














-1 









and the Hamiltonian vector field 



X = PdHi 



dV{qi,q2) dV{qi,q2) 



dqi 



We suppose that this vector field X is the bi-Hamiltonian vector field with respect to a 
certain Poisson bivector P', and some second integral of motion H2 

X = PdHi = P'dH2. 

The bivector P' has to be compatible with canonical bivector P, and it has to be the Poisson 
bivector. Hence, we have the following relations 



IP,P'1=0, [[P',P'J=0, and {Hi,H2} = {Hi,H2}' = 0. 



(3.7) 



Additionally we postulate that bivector P' is a Lie derivative of canonical bivector P along 
the vector field Y proportional to the Hamiltonian vector field X 



where multiplier p is equal to 



P' = CyP, Y = pX, 



p = apiqi + ^P2q2, 



(3.8) 



and so it coincides with N in (|1.5|) . 

Assumption (I3.8D ensures that P and P' are compatible bivectors, i.e., JP, P'] = 0, and 
that eigenvalues of the corresponding recursion operator M = P'P^^ are functions which 
depend only on Hi . In fact, we have 



det(A/'-AId4) = A2(^A + (a + |6)Hij , 
see (|3.6|) . Thus, we have constructed an irregular bi-Hamiltonian manifold M^. 
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Theorem 3.1. Bivector P' (|3.8|) satisfies the Jacobi condition 

|P',P'1=0 

c = o, V= p = /5p2<?2 



or 



where f and g are arbitrary smooth functions, a 7^ 0, /S and 7 are arbitrary numbers. 
Proof. In order to prove this theorem we have to substitute the second bivector 



qi 

(a+/3)/a 



(3.9) 

(3.10) 
(3.11) 



P' = CyP 



-locipiqi+l^p2qi -lapipi + ocqi 







dV 
dV 



dqi 



dV 
dV 



(3.12) 







dqi 

dV , dV 
-ocpi- hj6p2 5— 



into the Schouten bracket (|3.9|), and to solve the resulting system of partial differential 
equations. □ 



Theorem 3.2, Hamilton function Hi defined by dl.lD zfff/i potential V of the form (|3.10|) or (13.111) 

satisfies the following relations 

{Hi,p} = -(a + /3)Hi, (3.13) 



-(a + /3)Hi 



and 



(3.14) 
(3.15) 



X = PdHi = -(a + j6)"^P'dInHi. 

This therem is a direct consequence of the special form of P' and Hi . 

Inspired by the equations (|3.5|l for the Calogero-Moser super-tntegrable system we sup- 
pose that our systems are bi-Hamiltonian systems with respect to logarithm of the second 
integral of motion H2, i.e., 

X = PdHi = K2 ^P'd In H2. (3.16) 

This additional assumption is equivalent to the requirement that H2 and p satisfy supple- 
mentary relations 



{H2,p} = K2H2, and {H2,p}' = -(a + ^)k2 H1H2 . 



(3.17) 



The obtained relations (13.131) , (13.141) and (I3.17|) show that generators Hj, H2 and p form 
linear and quadratic Poisson algebras with respect to the brackets {■, ■} and •}', respec- 
tively. 
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Remark 3.3 Assumptions (|3.8|l and (|3.16|) yield the representation of the algebra l|1.2|) , but we can 
derive 113.161) /row algebraic relations (11.21) on/y (13.811 holds. Of course, form of (13.1611 a little 
strange because the logarithm of an integral of motion is itself integral of motion. 

The equation (|3.16ll is equivalent to the following system of the partial differential equa- 
tions 



K2H2 



2p2,2pi,-— — 



T 



aH2 aH2 aH2 dH2 



T 



(3.18) 



' dq2 ' dpi ' dp2 

with respect to unknown function H2 = H2{pi, P2, cii,cii), parameters a, j6, K2 and functions 
ficji) or ^((^j ^^'^^i) which appear in definition of potential V{q\,q2) given by (|3.10|) , or 
(13.111) . The explicit form of P' appearing in (13.181) is given by (I3.12|) . 

Remark 3.4 Sometimes equations (|3.18|) have different functionally independent solutions H2 as- 
sociated with the same potential V. In this case we obtain the so-called super-integrable system 
j|2¥ll25l/ , and our algebra (|1.2|) is some subalgebra of the complete polynomial algebra of integrals of 
motion. 

In order to get integral of motion H2 we additionally assume that it is polynomial in 
momenta, i.e., it has the form 

M 

^^2= E Shn{qi,qi)ViVi (3-19) 

k,m=0 

where gkmi^ir^i) are smooth functions. Substituting this expression into (|3.18|) we obtain 
partial differential equations for unknown functions gkf,i, f, g, and parameters a, ji. We can 
solve them for a fixed value of M. For instance, in the cubic case M = 3, for a = and 
K2 = 0, we obtain the following solution 

(Cl(?l + C2)"^ „ „ fCiqi + C2 

P = -qivi, ^ = - — — / H2 = Hi — ^ pi+p 

qi \ ci 

Of course, the Hamiltonian 

TT o , {ciqi + ci)'^ 

Hi = 2pip2 + — —, 

qi 

possesses the first integral 

Ciqi + C2 

h = — Pi - qivi, 

Ci 

which is linear in the momenta. However, this first integral I2 does not satisfy (|3.18|) in 
contrast with cubic integral H2. 

In the next step, if we look for a first integral of degree M = 4 in the momenta, then for 
a = 0, we reproduce previous solution and get one new 

p = -p2q2, y = —, f = {ciql + C2qi+C3y'^^^, 

qi 

associated with the fourth order in momenta integral of motion 



The obtained potential V is non-polynomial and non-homogeneous function. Below we 
will select from generic solutions certain particular solutions which are given by homoge- 
neous functions. 

4 Bi-Hamiltonian systems with second order integrals of motion 

Let us consider potential ^((Ji, qj) given by (I3.11II with 

7 = 0, and g(z) = . (4.1) 

That is 

and so, in this case the Hamiltonian has the form 

Hi=2pip2 + (?i " 4 (4.2) 

Taking this Hamiltonian we look for solutions H2 of equations (|3.18|) . We say that a solu- 
tion of these equations is trivial if the corresponding V[(]\,c\'i) does not depend on both 
variables. We do not distinguish solutions which are obtained one from the other by a 
permutation of pairs of variables and {fi,c\7). 

Substituting Hamilton function Hi of the form (|4.2|) into the equations (|3.18|) , and taking 
H2 of the form (I3.19|) with M = 2, one can get the following result. 

Theorem 4.1. if K2 = 1, a 7^ and j6 7^ equations (|3.18|) have only four solutions V{qi,q2) 
related to polynomial integrals of motion of degree two in the momenta 



1 


a — 




Hi 


= 2pip2- 


^q\qi 


h 






H2 = 





2 


a = 


-2^{d + l) 


Hi 


= 2pip2 - 


^ 4 


h 


= ^piiqipi - piqi) 




H2 = 




3 


a = 


2 


Hi 


= 2pip2 - 


Vqiqi 


h 






H2 = 




4 


a = 




Hi 


= 2pip2 - 




h 


= {piqi - pimf - 


2qi+^ 


H2 = 





Remark 4.2 In the case 2 and 3 H2 = h satisfies equations (|3.18|) with K2 = —cc, and K2 = —2a, 
respectively. Similarly, in case 4 H2 = l|*/Hi satisfies equations (|3.18ll with K2 = 2cc. 

For each case in the above proposition we suppose that h = Hi and I2 are the ac- 
tion variables, and co^, CO2 are the corresponding angle variables. Of course, this must be 
checked case by case. 
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In these variables the equations of motion have the form 

= 0, 12 = 0, cbi = l C02 = 0. 

Thus, C02 is a first integral functionally independent with Ii and I2 Il23| . The algebraic 
relations (I3.13II , (13.141) , and (|3.17|) allow us to express p in term of the action-angle variables. 
In fact, we have 



for case 1; 




= 0, 


P = 






for case 2: 


{PJ2} 


= a^2/ 


P = 




h j6)Jia;i - ttJ2(^2 + f (^1, h), 


for case 3: 


{PJ2} 


= loch 


P = 






for case 4: 


{PJ2} 


= 


P = 




V^)hcoi + F{h,l2). 



(4.3) 



Here angle variables cok are defined up to canonical transformations cok — > (jO]^ + f{lj^), for 
k = 1,2. 

Later we show that families of systems given by Theorem 14.11 are super-integrable for 
specific values of parameter d. That is, for these values of d the systems admit three 
functionally independent first integrals Hi, I2 and H3. Some examples of super-integrable 
systems with additional integral H3 of third, fourth and sixth order in momenta may be 
found in |l23l|2l. 

According to 11231 l24l 125 II , additional integral of motion H3 is a function on the action 
variables h, I2 and one angle variable co2- Below we prove that if we know additional 
polynomial integral of motion H3, then, using complete algebra of integrals of motion, we 
can get the angle variable CO2 algebraically. 

Some examples of super-integrable systems with additional integral H3 of third, fourth 
and sixth order in momenta may be found in 11231 124]| . 

Remark 4.3 Usually, see, e.g., [HI |2ll [251/ , in the theory of super-integrable systems we study 
polynomial Poisson algebra of integrals of motion. In this paper we add some extra generator p to 
this algebra. Let us recall that p plays the central role in our bi-Hamiltonian construction (|3.8|) . 

5 Case 1. Radial potential 

In this section we consider Hamiltonian system corresponding to case 1 in Theorem 14.11 
and in particular we look for values of d for which the system is super-integrable. 

We introduce new canonical variables (r, cp, p,-, Pq,) defined by the following equations 

(ji = r(cos ^ — i sin^) = re^'*', (^2 = f (cos ^ -|- i sin ^) = re"'', (5.1) 

and 

e'V / i \ e-'V / i \ 

Hamiltonian Hi and first integral /2 := — i/2 in new variables read 

Hi = ^ + ^ + V{r), V{r)=r^', J2 = p,. (5.3) 
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Thus after this transformation we obtain Hamiltonian of the classical problem of the 
motion in the field of the radial force written in polar coordinates. In the Cartesian coordi- 
nates x\= r cos (f, and xi = r sin q>, Hamiltonian Hi and the first integral /2 read 

= + r2 = x? + 4 /2 = yiX2 - y2Xi, (5.4) 

where j/i and \ji are the momenta conjugated with X\ and Xi. Thus we obtain a natural 
Hamiltonian system with the standard form of the kinetic energy and the radial potential. 
The problem of the existence of one more functionally independent first integral of this 
system is the question about maximal super-integrability. For natural Hamiltonian systems 
with homogeneous potentials possessing a proper Darboux point necessary conditions of 
the maximal super-integrability with first integrals which are meromorphic functions of 
coordinates and momenta were formulated in lfT2]| . These conditions were obtained from 
the analysis of invariants of the differential Galois group of variational equations. The 
considered radial potential possesses infinitely many proper Darboux points and results of 
IIT2II can be directly applied. In fact, in Section 3 of this paper Hamiltonian given by Hi in 
(|5.4|) was considered and the following result was proved. 

Theorem 5.1. The radial -potential (15.311 is super-integrable iff d = —1/2 or d = 1. 

Both of distinguished cases are indeed super-integrable with the additional first inte- 
grals 

H3 = (pi - P2) {piqi - piqi) - l^t^ , (5.5) 

for d — —1/2, and 

H3 = 2pl + ql (5.6) 

for d = 1, respectively. This result reminds us the classical Bertrand theorem |l2l which 
states that the only radial potentials V{r) = ocr^, for which all bounded orbits are periodic, 
are those with k = —1 and k = 2. The condition that all bounded orbits of a Hamiltonian 
system are periodic means that the system is degenerated, i.e., all its invariant tori are one 
dimensional. Such degeneration appears if the system is maximally super-integrable. 

6 Case 2 

In this section we consider Hamiltonian system corresponding to case 2 in Theorem 14.11 
In particular, we distinguish all values of d for which the system is super-integrable. We 
achieve this thanks to partial separation of variables and a direct integration of equations 
of motion. In IITtII this system is called quasi-separable and in fact we were able to find 
only one separation relation. 

6.1 Algebraic super-integrability 

One can easy observe that the syzygi exists between Hi and h 

h = -2p\qi + qiHi, (6.1) 
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but this is not yet a separation relation. We have to make the change of independent 
variable t — > t defined by dr/ dt = 1/ ^/qi, that transforms our Hamiltonian system into 

where ' denotes differentiation with respect to t. This transformation changes momenta 
and now we obtain one separation relation 

1 

It gives 

where is a constant of integration. Here we chosen square root with sign +, that corre- 
spond pi > 0. This integral is elementary and we obtain 



Hfy2((j2Hi- J2) = T + /5i. (6.3) 
In further calculations we choose j6i = 0, and then we obtain 

2/2 + H}t^ , ^ 

Next, we substitute this function into integral Hi, that is, into equation 

H = 1 dqid^2 ^ qj 
^ Iqx dr dr 

and obtain 

2^d7 + l2Hi)Vr= 
Then, we introduce new dependent variable v = v{x) defined by 

Making this substitution in equation (|6.5|) and dividing it by Hi we obtain 



T di; (2J2 + n\Ty ^ ^ 

V dr (2/2)'^^ 



(6.6) 



Finally, we make the following change of independent variable 

H2t2 

12 



and from (16.61) we obtain 



2x^-v+(l-xY = 0. 
ax 



(6.7) 



This is a non-homogeneous linear equation. Dividing it by {1 — xY and differentiating 
with respect to x we obtain homogeneous equation of the second order 



,d^v 1- (l-2d)xdv d 

x(l - x)—r H ^; '- — -v = 0. 

^ 2 dx 2 



(6.8) 



In this equation we recognise immediately Gauss hypergeometric differential equation 

x(l - x)w" +[c- {a + b + l)x]w' - abw = 0, (6.9) 

with parameters 



1 

2' 



1 

'=2- 



The general solution of equation (I6.8D is following 



(6.10) 



where (^2, j^s are constants of integration, and 2Fi{a,b,c;x) denotes the hypergeometric 
function with parameters a, b and c. Let us recall that 2^1 b, c; x) is a solution of the hy- 
pergeometric differential equation (16.911 which is holomorphic at the origin, see Appendix. 
Substituting (|6.10|) into non-homogeneous equation (|6.7p we obtain 



0. 



Hence, fis = 1, and solution of equation (|6.7|) is given by 

v{x) = f>i\fx^ 2F1 



(6.11) 



(6.12) 



Next, we return to the original dependent and independent variables and express t as a 
function of q2, see (16.41) . After this transformation the solution of (I6.5D takes the following 
form 



(6.13) 



From the obtained expression one can easily find the following first integral 



H3 = V2P2 



(6.14) 



We prove the following theorem. 
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Theorem 6.1. Hamiltonian system given by 



is super-integrable with algebraic additional first integral (|6.14l) iff d = p, or d = —{2p — l)/2, 
for a positive integer p. 

Proof. The first integral (I6.14I) is algebraic iff the hypergeometric function 2^1 appearing 
in its definition is algebraic. This implies that the hypergeometric equation (|6.9|) has an 
algebraic solution holomorphic at the origin. 

For equation (|6.9|) the differences of exponents at singularities are 

1 1 

A = 1 — c = -, v = c — a — b = d + l, ji = b — a = —d -\- -. (6.15) 

This is why, by Lemma lA.ll see Appendix, the Gauss hypergeometric equation (|6.9|l is 
reducible. In fact, for an arbitrary value of d we have 

A + /^ + i/ = 2, -A + p + v = 1, A-p + v = ld + \, A + p-v = -Id. (6.16) 

Thus, always at least one of this expressions is odd integer, and this proves our claim. 

As equation (I6.9D is reducible its one solution is algebraic. In fact Wi = ^/x is a solution 
of this equation. As we remarked, if first integral (|6.14|) is algebraic, then equation {6.9} has 
a solution W2 which is algebraic and holomorphic at the origin. Obviously W2 is linearly 
independent with Wi . In effect all solutions of equation \6.9\ are algebraic. 

Now we apply Theorem IA.3I in order to find values of d for which equation (I6.9D has 
only algebraic solutions. The first assumption of this theorem is that all exponents at all 
singularities are rational. For equation (|6.9|) exponents at x = 0, ;c = 1, and x = 00 are given 
by 

0,^}, {0,d + l}, 

Thus, if all solutions of equation (|6.9|) are algebraic, then d must be rational. 

The next assumption of the theorem is that among four numbers (I6.16D exactly two or 
four are odd integers. This condition implies that d = p/2, with p G Z. 

Now we have to check the presence of logarithms in local solutions of (16.911 around 
particular singularities. At first we assume that p is even i.e. p = 2q and d = q ^ Z. In 
this case the difference of exponents is integer only for singularity x = 1. For q ^ N, the 
difference of exponents ism = q + l — — q + 1. Hence, according to Appendix, the set 
{m) is (m) = {1,2, . . . ,q -\- 1}, see (|A.5|) . Now we have to check if 

^(m), for /,/G{l,2}, (6.17) 

with = q + 1. But we see that q + 1 + — q = 1 G (m). Thus for d = p/2 = q 
positive integer solutions of (16.911 contain no logarithms. For q negative integer difference 
of exponents is m = — {q + 1) = —q — 1, so (m) = {1,2, . . . , —q — 1}. Now we have to 
check all conditions (|6.17|l with 614 = 0. We obtain successively + — 1/2 = —1/2 ^ (ra), 
O-^O-cj = -q ^ {m), + 1/2 -1/2 = ^ (m), and + 1/2 = -q + 1/2 ^ (m). This 
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means that if is a negative integer, then logarithmic terms appear in solutions of (|6.9|) . 
Similar considerations for p = 2q — 1 show that logarithmic terms appear in local solutions 
of (|6.9|) around x = ooiffd = p/2 = {2q — l)/2, with positive integer q. 

When hypergeometric function 2^1 is an algebraic function of its argument, then H3 
defined in (16.1411 is an algebraic function. 

Finally, one can check by direct substitution that for d = hypergeometric function 
with the above parameters is equal to 1, and this ends our proof. □ 



6.2 Additional integrals polynomial in momenta 

Let us notice that in general the first integral given by (16.1411 for values of d specified in 
Theorem 16. II is not polynomial in the momenta. However, we can show the following. 

Lemma 6.2. If d = p E N, then 

H3 = H3,p := HaHf, (6.18) 
with H3 given by (I6.14D is a first integral of the system generated by Hamiltonian 

V 

H^Hi^2piP2' '^^ 



/qi 

Moreover, H3 is a polynomial in momenta and has degree 2p + 1. 

Proof. Clearly H3 is a first integral of H. In order to show that H3 is a polynomial in 
momenta we notice that 



H3 = H3H[ = ^ 



where 



h 



Wp(x):= 2fi(-^,-p,^,xj, Wp(0) = l, 
is a polynomial of degree p, see Appendix. Thus, 

is a polynomial in momenta of degree 2p + 2. Moreover we have 

V^Hr = ^ + ..., and l|Wp^^^^-^ 



where dots denote terms of degree at lest one in pi. This shows that B is divisible by pi. 
In effect H3 is polynomial in the momenta and has degree 2p + 1. □ 

Lemma 6.3. If d = —{2p — l)/2for a certain p G N, then 

H3 = H3,p:=H3J|/Hi, (6.19) 
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with H3 and I2 given by (|6.14|) and l|6.H) zs a first integral of system generated by Hamiltonian 

H = Hi = 2pip2 



Moreover H3 is a polynomial in momenta and has degree 2p — 1. 
Proof. First we show that for the assumed values of d we have 

2Fi(^-l,-d,^,x^ = 2Fi(-i l(2p-l),l,^) =(l-x)-(2p-3)/2Wp_i(x), 

where Wp_i(x) is a polynomial of degree p — 1, and Wp_i(0) = 1. 

In fact, for an arbitrary p, hypergeometric equation (|6.9|) has in a neighbourhood of 
X = 1 a local solution of the form w = (1 — x)^^'^P^^^^'^u{x) where u{x) is a function 
holomorphic at the origin and u(0) = 1. This is so, because — (2p — 3) /2 is an exponent at 
X = 1 for this equation. In order to find u{x) we make substitution iv = (1 — x)^^^P^^^^^u{x) 
in equation (|6.9|) . The obtained equation is again a hypergeometric equation. Its solution 
holomorphic at the origin is 



u{x) 



2F1 ( 1,1 - p,-,x 



(6.20) 



As p is a positive integer, u{x) is a polynomial of degree p — 1- 
Now, we have 



H3 -.^H^ll/Hi 



As in the proof of the previous lemma it is easy to notice that expression in the square 
bracket is a polynomial in momenta which is divisible by pi, and, as a result, H3 is a 
polynomial with respect to the momenta of degree 2p — 1. □ 

6.3 Examples 

Here we consider two examples of super-integrable systems given by Theorem 16.11 Let us 
take d = p = 2. In this case the Hamiltonian function has the form 



H = Hi = 2pip2 



Then, according Theorem 14. 1[ first integral h is following 

I2 = 2pi(p2r?2 -Piqi) 



4 



^1 

The second additional first integral H3 = H3 p which guarantees the super-integrability is 
given by (I6.I8II . For p = 2, it has, up to a multiplicative constant, the following form 



H3 = q^' -I6plql^^ + 3P2^2 + ^^PiPiqiipi + ^^qi) + 12pi{pi^ql + qiqi) 
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Moreover, there exists other first integral H4 of the form 

It is functionally independent with Hi and l2- Integrals H-[, I2, H3 and H4 are algebraically 
dependent as we have the following identity 

H4 + H3 + 32jf = 0. 

We note that 12^^"", Hj^^''^''\ H^^^" satisfy system (|3l8)) with normalisation K2 = 1. 
Examples with other values of p allow us to formulate the following. 

Conjecture 6.4. In the case 2 of Theorem \4.1\ if d = p ^ N, then there exists ctp G M, such that 
first integral H4 defined by the following relation 

H[+^H4 + Hip + apjf +^ = 0, 

is polynomial in the momenta of degree 2p. Moreover, integrals Hi, I2 and H4 are functionally 
independent. 

Now, we consider the element of family with d = —{2p — l)/2 for p = 3. Then Hamil- 
ton is 

Hi = H = 2pip2 + 



/qiq2 

and first integral I2 takes the form 

1 

k = 2pi{p2q2 - piqi) + ^ 3/2 - 

V^T2 

The explicit form of additional first integral H3 = H3 p defined by (|6.19|) is, up to a multi- 
plicative constant, equal to 

fr .„2.„ „ „ „ ^3 ^Pi{5plql-l2pip2qiq2 + 3plql) 5pi^i - p2?2 

H3 = 4pi(pi(?i - p2(?2j 3/2 ^ —3 ■ 

One can check that 

H rn« „^^3 I -33pf + 90pip2^ir?2 - 9pi^i 32 

H4 - 18pi(pi(?i - p2(?2j H ^ 3/2 ^-3' 

V^ql ^2 

is also a first integral of degree 2p — 2 = 4 with respect to the momenta, and that the 
following relation is satisfied 

32HI - 9H| - 7f H4 = 0. 

Moreover, /j^^''"^, Hg"^'''^'*-', H^^^'^ satisfy system (13.181) with normalisation K2 = 1. 
As in the previous case, additional examples justify the following conjecture. 

Conjecture 6.5. In the case 2 of Theorem 14.11 if d = — (2p — l)/2 ivith p G N, then there exists 
dip G M, such that first integral H^ defined by the following relation 

Hi + cipHlp - llHi = 0, 

is polynomial in the momenta of degree 2p — 2. Moreover, integrals Hi, I2 and H4 are functionally 
independent. 



17 



7 Case 3. Super-integrability without additive separation of vari- 
ables 

In case 3 of Theorem 14. 1 1 we have 

i^ = -rf^' ^ = * (Pi'?! " ^) ' H^^2p^p2 + qiqi (7-1) 
and one solution of the equations (I3.18D is equal to 

H2 = q'/^\ l2 = p\ + Kq{+\ ^ = J^- (7-2) 

As in the previous section our aim is to distinguish cases when the system is super- 
integrable. 

7.1 Algebraic super-integrability 

In order to find the second additional first integral H3 we proceed as follows. From the 
integral I2 (I7.2D we directly obtain the first quadrature 

ft+,=-r^L=. (7.3) 

After calculation of this quadrature we have to substitute solution q2 = qiihj^i) into the 
Hamiltonian 

_ldqi dq2 A 

and solve the obtained quadrature with respect to q\. However to find q2 = qiit) we have 
to invert explicitly quadrature (|7.3|) , and it is unclear how to do this. By this reason we 
proceed in a different way. Using function p and algebraic relations (14.31) we can get the 
second angle variable 

_ p d-1 k p2 dx 
"^^'2^^ JTT 222 J 2yj2-?cx^+i ■ 



In generic case 002 is the multi-valued function on the whole phase space 

We want to distinguish values of d for that system in case 3 from Theorem [41] governed 
by Hamiltonian Hi given by Il7.ll) is super-integrable. Since Hi = Ii, thus CO2 is a first 
integral. So, the following function 



H3 = 4a(rf + 1)720)2 = 2{d + l)p + {d- l)aJi J 



12 dx 



yJl2-KX^+^ 



is also a first integral which we are going to analyse in details. 

In the integral presented in H3 we make the following substitution 

y = ^^x'^\ 
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which gives 

1 _ d 

dx _ 1 ih\''+'^ fy ^dy 



Let us denote 

d 



f^-^ = {d + l)yMy), (7.5) 

where v{y) is a certain function of y. Differentiating both sides of this equation with 
respect to y, after obvious simplifications we obtain the following linear non-homogeneous 
equation for the unknown function v{y) 

v+{d + l)yv'. (7.6) 



Multiplying it by y^l — y and differentiating once more we obtain 

y^l _ yy ^y") - -i-. = 0. (7.7) 

\d^\ Id + l'^J 2d + 2 ' 

This is the Gauss hypergeometric equation with parameters 

1,1 ,1 

a = -, b = - — -, c = 1 + 



2' d + 1' d + 1' 

General solution of (I7.7D has the form 



Ky) = fcy-"<^«' + ft .ft (1 1 + . (7.8) 

Substituting it into equation (I7.6D we find out that it is its solution iff /33 = 1. Without loss 
of the generality we can assume that f^2 = 0, and then in variable x considered integral is 

dx X /I 1 1 KX'^+l 

2F1 ^,^—r^ + 



^l^-KX^+^ ^/h \2'd + l' d + r h 
Thus our first integral takes the final form 

This first integral is functionally independent with Ji and l2- In general hypergeometric 
functions are not single-valued with sinularities of non-algebraic character. Thus, in order 
to obtain first integrals with "good" behaviour, we will restrict to cases when function 2^1 
is algebraic. In other words, we are looking for such values of d that equation (|7.7|l has a 
solution that is algebraic. At first let us note that by Lemma lA.ll equation (I7.7D is reducible 
for any value of d. Indeed, the differences of exponents at singularities of equation ^7.7^ are 

A = l-c=-^, v = c-«-fc = l, ^ = b-« = ^-l, (7.10) 
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and among numbers 



\ + H + v = 0, -A + }i + v = j-^, + ^ = A + f/-i/ = -l, (7.11) 

we have one even and one odd number. This means that we can use Theorem I A.3I in order 
to find values of d for which equation (|7.7|) has only algebraic solutions. 

Theorem 7.1. Hamiltonian system given by 

is super-integrable with algebraic additional first integral (|7.9|) iffd takes the form 

d='-^, or d = \±^, (7.12) 
p 1 — 2p 

for p eN. 

Proof. We apply directly Theorem IA.3I Exponents at singularities x = 0, x = 1 and x = oo 
for equation (|7.7|) are 

11 f 11 fl 1 



d + lj' I '2J ' \2' d + 1 

respectively. The first condition of Theorem IA.3I implies that d has to be rational. The 
second condition can be fulfilled only when two among number l|7.1Hl are odd, and this 
implies that 

2 2-r 

-y—r = r, so d= , (7.13) 

a + 1 r 

for a certain r G Z. 

Now we have to check the presence of logarithms in local solutions of (|7.7|) around 
particular singularities. In this aim we apply Lemma IA.2I For d given by (17.131) exponents 
are following 



»-^}' {M} 



We see that differences of exponents can be integer positive only for x = or x = oo 
depending on the parity of r as well as the positive or negative sign of r. An analysis 
similar to that done in the proof of Theorem 16.11 shows that logarithmic terms appear in 
local solutions iff either r is positive and odd, or r is negative and even. 

If all solutions of the considered equation are algebraic then, by Theorem IA.3I none 
of local solution has a logarithmic term. Thus, by the above considerations, either r is a 
positive and even integer, or it is a negative and odd integer. Substituting r = 2p as well as 
r = —{2p — 1) with p G N into (I7.13|) we obtain (|7.12l) , and this ends the proof. □ 
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7.2 Additional first integrals polynomial in the momenta 

If assumptions of Theorem 17.11 are fulfilled, then first integral given by II7.9I) is algebraic 
but it is not polynomial with respect to the momenta. However we can find such integrals 
using the following. 

Lemma 7.2. If 

1-p 

d = for a certain p G N, 

P 

then H3 given by (I7.9D can be uniquely written in the following form 

H3 := H3,p + /?2/3,p, (7.14) 

where Hs^p and /a^p are first integrals polynomial in the momenta of the system with Hamiltonian 
H = Hi defined by (I7.1I) . Moreover, Hs^p has degree 2p — 1 with respect to the moment and is 
functionally independent together with H and h- 

Proof. For the specified form of d, we have 

F{x):= 1 + = 2F1 [1,PA + P,xy (7.15) 

We show that 

F{x) = x-P (C + yr^Wp^i (x)) , (7.16) 
where Wp_i (x) = 2F1 {■^A — pA — p, x) is a polynomial of degree p — 1, and 

c-2F.[rP^^^pA)-^Yir^-j^^: 

In fact, V := F{x) is a solution of hypergeometric equation (|7.7|) holomorphic at the origin. 
Making substitution (17.161) we find that Wp_i (x) coincides with 2^1 {\A — p,'^ — p,x). As 
p is a positive integer Wp_i(x) is a polynomial of degree p — 1. 
In the considered case we have 

Thus we can rewrite H3 given by \7.9\ in the following form 



i.e., it has the form (|7.14|) with 



H3,p := a^(pii?i - 2pp2^2) + a^-^pih^ Vp_i (^^^ j , (7.18) 
and 

hp ■■= «C^Iir'- (7.19) 
Clearly, these are polynomials with respect to momenta, and Hs^p has degree 2p — 1. □ 
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Lemma 7.3. If 



d = for a certain p G N, (7.20) 



then 

tP 



H3,p := Hsq, (7.21) 
where H3 given by (I7.9D is a first integral polynomial in the momenta of degree 2p + 1. 

Proof. We show that if d is given by (|7.20|l , then 



22) 



where Rp_i(y) is a polynomial of degree p — 1. 

Recall that hypergeometric function given above is a solution of equation (I7.7D . If we 
make substitution 

v{y) = v^r^u(y) (7.23) 
in equation (|7.7|) , then we obtain again the Gauss hypergeometric equation 

y(l - y)u" +(^^-p + {p- 3)y^ w' + (p - !)« = 0, (7.24) 

with the parameters 

a = l-p, b = l, c = ^(3-2p). (7.25) 

For p G N, parameter a is non-positive integer, thus function u{y) := 2F1 (,i,b,c,y) is a 
polynomial of degree p — 1. Hence, we have 

/II 3 \ 

v{y) = 2F1 ( 2'2~^'2 "'^'^J " -]/"(!/) = 

yr^2fi (l-p, 1,^(3 -2p),y) . (7.26) 

It follows that Rp_i(y) is equal to 2F1 (l — p,l,\{3 — 2p),y), and so it is a polynomial of 
degree p — 1, as we claimed. 

Now we have to show that Hs^p given by (|7.21|) is a polynomial with respect to the 
momenta. But we can write 



H3,p = cell {l{d + l)^?ipi - q2p2) + «(d - m^-^J 1 - ( 



(2(d + l)^jipi - q2p2) + «(d - l)hqiVT2\l 1 - ^^Sp_i, (7.27) 



where 
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J2 



is a polynomial of degree not higher than 2{p — 1) with respect to the momenta. Now, it is 
enough to notice that according to (I7.3II and (I7.2D , we have 



in order to see that both terms in (|7.27|) are polynomial with respect to the momenta. 
Moreover, it is easy to see that the degree of H3 p with respect to the momenta is 2p + 1. □ 



7.3 Examples 

As example we take the family with d = (1 — p)/p, and we choose p = 5. In this case we 
have 

Hi=2pip2 + 475' h = pl+5cjY'. 

Up to a multiplicative constant, the first integral H3 := p defined in Lemma 17.21 by 
equation (|7.14|) has the following form 

H3 =^?2 ^/^[128p?p2'?2^^ + 64p^(Pi'?i + 35p2^?2) + 560p\qY^{2p,q^ + 25p2'?2) 
+ 7000p]ql^^{piqi+5p2q2)+'^375ql^^{4:p^qi+5p2q2)]- 

There exists also other first integral of degree 2p — 2 = 8 in the momenta 

Hi =S96p\p2{piqi - piqi) + llAq^^^^ p\{lp\ql + 66pip2qiq2 - 65plql) 
+ 16q2^^^pl{^76plql + 5215pip2'?i'?2 - SOlSpjql) 
+ 5q2^^^{9072plql + 32920pip2^?i(?2 - 30625p|(?|) + 102400(??(?2 

and the following relation holds 

3125H4 + iff - mGH^ll = 0. 

Moreover, Jj^^^^^'*^ Hy^'^'^\ h]^^'^'^' satisfy system (I3.18|) with normalisation K2 = 1. Other 
examples justify the following. 

Conjecture 7.4. In the case 3 of Theorem \4.1\ if d = (1 — p)/p for a certain p G N, then there 
exists a.p G M, swc/z that first integral H4 gfuen fcy 

is polynomial in the momenta of degree 2p — 2. Moreover, integrals Hi, I2 and H4 are functionally 
independent. 

Now, let d = (1 + 2p)/ (1 — 2p), and we consider the case with p = 2, i.e., d = —5/3. In 
this case we have 

Hi=2pip2 + 473' and 72 = P? - -I73 • (7.29) 

'?2 ^r2 
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The additional first integral H3 := H^^p defined by (|7.21|l of degree 2/3 + 1 with respect to 
the momenta is, up to a multiplicative constant, of the following form 

ffs - ¥Apm - pm) + "^'^""%7 ""''^ - ^'^P-'" + ^"'"^^ . (7.30) 

E 4 

We have also first integral H4 of degree 2p with respect to the momenta given by 

rr_9.^„ r, n ^2 (npi^l-27p2'?2)(pi?l-p2q2) 16^? , . 

2^J2 Cj2 

and the following identity holds 

864H? - + 16H4jf = 0. 

Moreover, l2^^^'^'^\ H^^''^^'^\ H^^''^^"^ satisfy system (I3.18I) with normalisation K2 = 1. 
Other examples with different values of p allow us to formulate the following. 

Conjecture 7.5. In the case 3 of Theorem \4.1\ if d = (1 — 2p)/ (1 + 2p) for a certain p G N, then 
there exists Up G M, smc/z that first integral H4 defined by 

H^I^^^ + Hf^p + ftpHi^ = 

is polynomial in the momenta of degree 2p. Moreover, integrals Hi, I2 and H4 are functionally 
independent. 

Let us determine the complete algebra generated by first integrals (I7.29D , (I7.30D , (I7.31I) 
and p which in the case d = —5/3 has the following form 

1 

P = - 2(^1^1 +3'?2P2). 

It is easy to check that first integral given by 

H3 = {h,H4} , 

satisfies equations (|3.18|l with ^2 = 2. The complete algebra of integrals and p is given by 

{h,p}=2h, Jl2,p} = l2, {H4,p} = H4, {H3,p}=2H3, {I2,H4}=H3, 

{Ji, ^2} = {h, H3} = { Ji, H4} = 0, {I2, H3} = 8J|, {H3, H4} = 2472H4. 
Using this algebra we easily obtain the desired angle variable 



CO2 



up to canonical transformation CO2 — > (^2+ f{h)- Here F is arbitrary function of the con- 
stant of motion, chosen in such a way that relation {001,002} = is satisfied. 
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8 Case 4 

In case 4 of Theorem 14. 1 1 we have 



H^=2p^p2 + q^'-y2- (8-1) 

The system generated by this Hamiltonian is integrable for an arbitrary value of parameter 
d with the following first integral 

l2 = (pi^i-p2^2)'-24^. (8.2) 

This system admits separation of variables and we show that for rational values of d it is 
super-integrable. 

8.1 Super-integrability and additive separation of variables 

Let us introduce the following canonical coordinates (f, <p) and momenta {pr,P(p) defined 
by 

qi = r(cosh cp + sinh cp) = re^ , q2 = r(cosh cp — sinh cp) = re ^ , (8.3) 



and 

~ " 2 V ' r J ' 2 r 

In new variables Hamiltonian Hi and first integral = li/2 take the forms 



P. = ^(^ + ^). P2 = ?fp.-^). (8.4) 



2 2r2 



Hi = ^-^ + ^—: /2 = ^pJ-e-^^+^^^. (8.5) 



New Hamilton's equations read 

-p2 +2e-2('^+i)«' 



r = Vr, Vr ~ ^3 



2(d + l)e-2('^+i)'P 

Let us note that 



(8.6) 



Hi = ^p?-^/2. (8.7) 

In order to perform the explicit integration we introduce as in fH a new independent 
variable t such that dr/dt = Then we find that 

where prime denotes the differentiation with respect to t. In effect we have 
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i.e., we effectively separated variables and we can make two independent quadratures 



dr 



T+Ci, 



d(p 



2(/2 + e-2(rf+i)9') 



T + C2. 



The explicit forms of these elementary integrals are following 



1 



In- 



and 



1 



(d + l)v/2]2 

Using well-known formula 



arcsinh 



{d+l)(p 



T + Ci, 



T + C2. 



arcsinh z = ln[z + vz^ + l], 
one can rewrite the last integral as 

1 



In 



From §^ and dSlOl) we deduce that 

7= (d + l)y2/2(C2-Ci) 
= ln 



T + C2. 



(d + l)ln 



In 



2(v^+^/2 + Hir2 
7/2 + Je^id+^)vh + 1 ) ( v//2 + V/2 + Hir2' 



is a first integral of the system. We use relations 



h + Hir^ 



2 2 



p2(rf+l) 



^2 + 1 



2 ' 2 

to introduce explicitly the momenta and then I takes the form 

-2'i/2 



2(rf+l)<p 2 



; = In 



rd+1 



W2j2-p^)W2j2 + rpr 



\d+l 



We can take as the first integral its exponent, more precisely 

H3 = 2-'*/2exp(J) = -l_e('^+i)^'(^_p^)(y]2 + rp,,)'+'. 
When we come back to original variables, then this first integral takes the form 

H3 = -j^ipicji - piqi + vh) ipiqi + piqi + \/h\ 

^2 



(8.8) 



(8.9) 



(8.10) 



(8.11) 



(8.12) 
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It is easy to notice that 

1 / ^\'*+^ 

H4 = -^{fiqi - piqi - vh) ipiqi + piqi -Vhj 

is also a first integral. For rational d these functions are algebraic. Thus we can recapitulate 
these considerations by the following theorem. 

Theorem 8.1. Hamiltonian system given by (18. ID is super-integrable with algebraic additional first 
integral iff d is a rational number. 

Using H3 and H4 we can get first integrals which are polynomial with respect to the 
momenta. For example, if d is positive integer, then we take the following first integrals 

_ 1 

H3 = —r={H3 — H4), H4 = H3 + H4. 



In the sum H3 + H4 terms containing odd powers of i/ia disappear, and, as a result, H4 is a 
polynomial with respect to the momenta of degree d + 2. The difference H3 — H4 contains 
only odd powers of 1/I2 , thus it is divisible by 1/I2, and H3 is polynomial in the momenta 
of degree d + 1. 

If d is a positive rational number of the form d = di/ d2, then we put 

F4 = = -j^{p2q2 - piqi - VhY' [piqi + P2^?2 - ^/h) ^ , 
^2 



and then 



_ 1 

H3 = ^=(f3 — f4), H4 = f 3 + F4, 



are polynomial first integrals of degrees d\ + 2.di — 1 and d\ + 2.di, respectively. 

In the case when d > — 1 is rational negative, i.e. d = —d\ldi, and ^2 > d\ where 
dx,d2 G N we take 

H3 = 4^ (Hg^^ - 2 ), H4 = Hg^^ + H{\ 

as polynomial in the momenta first integrals of degrees 2d2 — di — 1 and 2^2 — di, respec- 
tively. 

It is not clear if there exist a polynomial in momenta first integral functionally indepen- 
dent with Hi and h for rational d smaller than —1. 

Let us note that potentials from case 4 written in Cartesian coordinates as 

(x - iyf+^ 

were considered in the recent paper [51. In this paper authors proved the super-integrability 
of such potentials for k rational for classical and quantum systems. But the explicit form of 
the additional first integral was not announced. 
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8.2 Examples 

We showed that for rational d the system is super-integrable. We have the second addi- 
tional first integral H3 which can be chosen as a polynomial with respect to the momenta. 
Moreover, there exists another additional polynomial integral of motion H4 which can be 
normalised in such a way that 

Ha = {^2/^3} , 

satisfying the equations (I3.18D as well. For K2 = 1 adding polynomial integrals to the 
generators p, h, h we obtain the complete algebra of integrals of motion and function p 
defined by the following brackets 

{h,p} = -2(xh, {h,l2} = 0, {h,Hk} = 0, {p,l2}=0, 

{I2,H3} = Ha, {H3A'P} = H3A' {I2,Ha} = 2k^I2H3, {Ha,H3} = k^HI. 

(8.13) 

Here parameter k depends on a and d. Let us present some particular examples with cubic 
polynomial integrals 

. iora = --,K = 18,V=^: 
3 

H3 = 2pl - 2tj2P2^jr^ - Pi^i^' 

~ bet' 

H4 = I2p{{p2qi - piqi) - (^qi^iplql + pipiqiqi + ^piqi) + 

qi 

36Hf -H| + 36H^72 = 0, 
1 

• for Di = —1, K = 2,V 



ql'^qr 



H3 = 2p2(m2-m0^-^q^, 

qi qi 



q^ql" q\"qr' 

fourth order integrals 

• for a = -\,K = 32,V =% 
4 q\ 

rj _ 4 ^qlpl qipipi p\ q\ 

' 2q\ ql 2q\^m' 

u ^^'^(^ a 4(Pi'?i + PiP2?i^2 + Pipit?itj^ + 5p^tj^^ 2qj{3piqi + 5p2q2) 

W4 = bp2[p2q2 - piqi) 4 \ 8 ' 

qt qi 

8Hf -Hl + (AHlh = 0, 
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for a = -1,K = 2,V 



1-13 - Plipiqi - piqi) +7j-7=z ^ 3/4 1/4 



H4 = -2p2ipiqi - piqi) 



iql'W' 

4 - 5p2q2){piqi - piqif ^piqi - 5p2t?2 



and fifth order additional integral of motion 
• for a = --, K = 50,V =^ 

5 

H3 = -p2 (ipl - f2pi - - ^ - ^'^^P^Pl + '?2Pl('?2 - ^Plql) 



2 V qlJ \ ql J q\ q{ q{ 

£5 7n«5/„ , „ „ ^ 10(p1?1 + P2?^t?2 + p?pi'??qi + Piplqiql + 6p|^|) 

j-;4 = i{jp2\piqi — piqi) 5 

qi 

5ql{^p\q\ + 7pip2qiq2 + 9p2t?2) _ 5^ 

q\° q? ' 

25Hl - H| + 100H|J2 = 0. 

According to 11231 125II H3 has to be function on the both action variables Ji 2 and the 
single angle variable C02- Using the complete algebra (|8.13|l we can easily get 



Here CV2 is defined up to canonical transformation LO2 — > c<^2 +/(^2) and parameter c is 
calculated from the bracket {l2/<^2} = 1- For instance, for a = —1/3 and d = —4 we have 
c = -l/16. 

On the other hand, if H3 is the known polynomial solution of (I3.18I) , which satisfies the 
algebra (|8.13|) , then we can directly calculate cv2 



1 / Hsh + ^Hiq-AchlY^ 



^2 = -^ In I '—-^^ I (8.14) 



without separation of variables. 

9 Bi-Hamiltonian systems with higher order integrals of motion 

Let us recall that in Section |4] we investigated systems given by multi-parameter family of 
Hamiltonian functions given by 

/i(ri+l)+a . 

H = 2pip2 + ^i " qi (9.1) 
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see (|4.2|) . For these systems with arbitrary values of parameters a, and d, the correspond- 
ing bivector P' given by (13.121) is the Poisson bivector compatible, by Proposition I3.1[ with 
canonical bivector P. In previous sections we distinguished infinitely many integrable cases 
with additional first integrals which are polynomial in the momenta of the second degree. 
Among them infinitely many systems are super-integrable. 

In this section we continue our integrability analysis and we look for cases when the 
additional first integral is a polynomial in momenta of degree greater than two. To this end 
we solve equations (|3.18|) with respect to H2 of the form (|3.19|) . The existence of solutions 
of the obtained equations depends on values of parameters oc, /S, d and K2. 

Assuming that the additional first integral is of degree 4 in the momenta and setting 
K2 = 1 we found the following solutions 

P = -^{qiPi+^piqi), V = qlqf/^, 

H2 = 4p4 - 10 {3plql - 30pip2qiq2 + Ibplql) q^^'^ + 215q\qf^ , 

9 = - -^{^p\q\-piqi) , y = qlqL 

H2 = \6pl{piqi - p2q2) + ^Piqiq\{piqi - 2p2'?2) + q\ (pi - ^hi) ' 

p = -^{piqi+^^piqi), y = qlq2^^^, 

H2 = Iplipiqi - piqi) - q2^'^ {^^phl " ^^PiPiqiqi + 64pi^?i) + GAqlq^^/^, 

p = -piqi-^piqi, ^^^i^^^il^*' 

H2 = ^Piipiqi - Piqif - 2qf^^qf^^ {3plq\ - I2pip2qiqi + plql) + 9q^^'^qf'^ . 
Assuming that H2 is of sixth degree with respect to the momenta we obtained other variety 
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of solutions: 

P = -l ^Viqi + 7P2^?2) , V = q\q2^'"\ 

H2 = 4pf - lOSSqlql^^ + 343(?2 - 7p2^2)(pi^?i - P2^?2) + Upl{-Ap^q^ + 7p2q2)q2 

p = -Piqi + '2.piq2, V ^ q-^'^q^^'^ , 

2/3 „ „ ^ ^2 / 1/6 T„ ^2/3/„ ^ ^ ^ ^^ , ^-1/3^-1/3 



H2 = 2^?^ - p2q2) [q'2: - 2piq{"{piqi - piqi)) + q^ '"q2 

p = -5(^1^1 +p2'?2), ■^^ = ^?^v^^^ 

R - « „6 „ „ „5 ^_i^.Hv\ , 3^ipiP2 , pi ^1^' 



2^^/2 2<j5 4^?4 4^3 s^l5/2 ■ 

Here we present one super-integrable system for which equation p.l8|l has two func- 
tionally independent solutions H2 and H3: 

p = -piqi - Piqi , V = qi^'^q2^'^ ■ 



In this case, for K2 = 1, there is fourth order polynomial integral of motion 

o _ viiviqi-nqif i3pf^f-44pip2?i?2 + 4pi^i , ,-1/3,-8/3 

and for K2 = 2 sixth order polynomial first integral 

H3 = 4p?(pi(^i-p2i?2)^ + (?r^^^'?2^^Ml0pl'?l-16piP2'?l^j2 + pi'?i) 
- 4:q^^^^q:^'^^^pi{piqi - p2q2){'2.p\ql - 6pip2qiq2 + plqj 



qiq\ 



In this case the complete algebra of integrals and p is given by 

{p,Hi} = (a + j6)Hi, {p,H2} = -H2, {p,H^} = -2H^, {p,H^} = -3Hi, 

3 27 
{H2,H3} = H4, {H2,H4} = - H3, {H3,H^} = -^H-^ ■ 

Remark 9.1 According to our knowledge it is the first example of super-integrable system with first 
integrals of degree two, four and six with respect to the momenta. In all examples known to us 
super-integrable systems have at least tivo additional first integrals of degree two with respect to the 
momenta |i [H IHl El ESP. 
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9.1 Additive separation of variables without super-integrability 

Let us consider one of the bi-Hamiltonian systems Usted above 



p = -^{qiPi+5p2q2) , K2 = l, 



H2 = - 10 {3v\q\ - 30pip2qiq2 + 25plq^) q^^^^ + llSqfqf^^ . 
In this case second Poisson bivector is equal to 



qipi-5q2p2 PiPi 



P' 



CpxP = 2 



2qr 



15£ 
2qr 



^q^ 

^^{3qipi + 25q2p2) 



lOq 



and the corresponding recursion operator J\f = P'P Ms degenerated 

det(Ar-AId) = a2 I^A-^Hi^ . 

According to p6], we have to look for another solution P' of the equations (|3.1|) associated 
with the non-trivial Darboux-Nijenhuis variables. The existence of such solution converts 
into a regular bi-Hamiltonian manifold. 

It is easily to find three different second order polynomial solutions of equations (|3.1|) . 
Only one of them allows us to get variables of separation. Namely, if 



and 



= -5(72 [i^PiPiql^' + qi)VlO - ?>piqi'^ql - I5p2qi'^qi 
are coordinates of the vector field Z = Yj Z^d^, then bivector 

P' = CzP 

is the solution of (13. ID for the given integrable system. For this Poisson bivector P' the 
control matrix F 

2 

P'dH; = P ^ Fij dHj, i = 1, 2, (9.2) 

7=1 

is a non-degenerate matrix of the form 

Wql/^{-15q^+ql^^VWpi) \ 
-500qY^{2q2pl + I5q\'^q\ - 3ql^^VT0piqi + 5ql^^VWp2) 



32 



The eigenvalues ui 2 of F 



A{A) = (A-Ui)(A-W2) = 

= - Wqy^{ql^^VWpi - 15qi)\ + 250ql^^(l5ql + VWql^^{5p2q2 - ?>qipi) +2v\q\'^) , 



coincide with eigenvalues of the corresponding recursion operator J\f = P' ■ P^^ — CzP ■ 
P^^. Thus Mi^2 are the variables of separation and the so-called Darboux-Nijenhuis coordi- 
nates. The conjugated momenta i^i 2 are equal to 

The separated equations look like 

- 2 {av^f Uk + H2 + liikHi = , k = l,2, a = IQ^'^ . (9.3) 

Thus, in order to get solutions 162) and Uk{t, jSi, 1^2) we have to solve the following 

Abel-Jacobi equations 



(t73t;3 + Hi)2 ^7 ((73^73 + Hi )2' 
/■"i du /''^2 dv 

^ J 2(c73y3-Hi) 2(t73z;3-Hi)' ^^'^^ 

There are no addition theorem for the above quadratures and, therefore, we can suppose, 
that there are no additional polynomial integrals of motion 1241 l25ll . In fact, the second 
angle variable 

1 ^ dv 1 2 

2^J aH^-Hi UaH^/^fi^ ^ 
- In [a'vl + av,H\'' + hI'') - 2V3arctan (^^^^|^^ 

cannot be rewritten as a function on polynomials Hi, H2 and H3 as in the super-integrable 
case (|8.14|) . 

In the similar manner we can get solutions of the equations of motion for other systems 
with fourth order integral of motion. Construction of the variables of separation for the 
systems with sixth order integrals of motion is an open problem. 
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A Appendix. Hypergeometric equation 

Hypergeometric equation 

z(l - z)w" + [c- {a + b + l)z]w' - abw = 0, (A.l) 

is a special case of the Riemann P equation. Thus it possesses three regular singularities at 
z = 0, z = 1 and z = oo, with exponents and their differences of the following form 



z = 0, 


pi 


= l-c. 


P2 


= 0, 


A = 1 - c. 


z = l, 


0-1 


= c — a — b, 




= 0, 


V = c — a — b, 


Z = 00, 




= h, 


Tl 


= a, 


ji = b — a. 



Exponents satisfy the Fuchs relation, i.e. their sum is equal to one. 

As it is well known one of solutions of equation jA.l} holomorphic at the origin is the 
hypergeometric function given by the following series 

,Fi(«,fc,c,z):=f;i^z^ (A.2) 

Here (x);^ = x{x + 1) • • ■ [x + k) is the Pochhammer symbol and {x)o = 1. The above series 
is a polynomial iff fl or is a non-positive integer, see, e.g., IITSl lTll. 

The following lemma gives the necessary and sufficient condition for (IA.1|) to be re- 
ducible. It is a classical, well known fact, see, e.g., llTll . 

Lemma A.l. Equation (|A.1|) is reducible if and only there exist i, j, k G {1,2}, such that 

Pi + aj + TkeZ. (A.3) 
Equivalently, equation HA.!)) is reducible if and only if at least one number among 

A + fi + v, —A + ^ + A — ]A- + v, k + n — V, (A.4) 

is an odd integer. 

From the above lemma it follows that if equation {A.l} is reducible, then we can always 
renumber exponents in such a way that 

pi + ai+Tie -No, 

where No denotes the set of non-negative integers. But then, from the Fuchs relation, we 
also have 

,02 + (7-2 + T2 e N. 

Hence, if dA.lD is reducible, we assume from now on that the exponents are numbered in 
this way. 

If the difference of exponents at a singular point is an integer, then it can happen 
that a local solution around this singularity contains a logarithm. Such a singularity is 
called logarithmic. For equation (|A.H) , it is enough to know the exponents to decide which 
singularity is logarithmic. To formulate the next lemma which gives the necessary and 
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sufficient conditions for a singularity of (|A.1|) to be logarithmic we introduce the following 
notation. For a non-negative integer m G No we define 

/ \ J® if m = 0, 

I{1, otherwise. 

For s G {0, 1, 00} let 65,1 and 65,2 denote exponents of equation dA.ll) , ordered in such a way 
that Re e,,! > Re £5,2- With the above notation we have the following. 

Lemma A.2. Let r G {0, 1, 00}. Then r is a logarithmic singularity of equation \A.l) if and only if 
m := e^,! — ^r,2 £ No, and 

er,i+es,i + et,j ^ {m), for i,j G {1,2}, (A.6) 
w/iere r, s, t are pairwise different elements 0/ {0, 1, 00}. 

For more details, see Lemma 4.7 and its proof on pp. 91-93 in Q. 

For considerations of this paper we need characterisation when reducible hypergeomet- 
ric equation has algebraic solutions. The answer is contained in the following theorem. 

Theorem A.3. Assume that hyipergeometric equation \A.l) is reducible. Then all its solutions are 
algebraic iff 

1. all exponents are rational, and 

2. exactly two or four of A + p + v, —A + p + v, A — p + v, A + p — vare odd integers, and 

3. no one of singularities is logarithmic. 
For proof and details see ll9l[7l[T8ll. 

References 

[1] M. Audin. Les systhnes hamiltoniens et leur integrabilite. Cours Specialises 8, Collection 
SMF. SMF et EDP Sciences, Paris, 2001. 

[2] J. Bertrand, Theoreme relatif au mouvement d'un point attire vers un centre fixe. C. 
R. Acad. Sci. Paris v. LXXVII, pp. 849-853, 1873. 

[3] O. I. Bogoyavlenskij. Algebraic properties of master symmetries and their applications. 
Dokl. Akad. Nauk, v. 360(4), pp. 445-447, 1998. 

[4] A. V. Borisov, A. A. Kilin, and I. S. Mamaev. Superintegrable Systems on a Sphere 
with Integral of Higher Degree. Regul. Chaotic Dyn., v.l4 (6), pp. 615-620, 2009. 

[5] C. Daskaloyannis, K. Ypsilantis, Unified treatment and classification of superintegrable 
systems with integrals quadratic in momenta on a two-dimensional manifold, /. Math. 



Phys., V.A7, 2006, p. 042904, 38 pages, arXiv:math-ph/ 0412055 y3 



35 



[6] Yu.A. Grigoryev, A.V. Tsiganov, Symbolic software for separation of variables in the 
Hamilton-Jacobi equation for the L-systems, Regul. Chaotic Dyn., v. 10(4), pp. 413^22, 
2005. 

[7] K. Iwasaki, H. Kimura, S. Shimomura, and M. Yoshida, From Gauss to Painleve, A 
modern theory of special functions, Aspects of Mathematics, E16. Friedr. Vieweg & Sohn, 
Braunschweig, 1991. 

[8] E. G. Kalnins, W. Miller Jr., G. S. Pogosyan, Superintegrability and higher order con- 
stants for classical and quantum systems, arXiv:0912.2278vl [math-ph], 24 pages. 

[9] T. Kimura, On Riemann's equations which are solvable by quadratures. Funkcial. 
Ekvac, V 12, pp. 269-281, 1969/1970. 

[10] A. J. Maciejewski, M. Przybylska, All meromorphically integrable 2D Hamiltonian 
systems with homogeneous potential of degree 3, Phys. Lett. A, v. 327(5-6), pp. 461- 
473, 2004. 

[11] A. J.Maciejewski, M. Przybylska, Darboux points and integrability of Hamiltonian 
systems with homogeneous polynomial potential, /. Math. Phys., v. 46(6), p. 062901, 33 
pages, 2005. 

[12] A. J. Maciejewski, M. Przybylska, H. Yoshida, Necessary conditions for super- 
integrability of Hamiltonian systems. Phys. Lett. A, v. 372(34), pp. 5581-5587, 2008. 

[13] A. J. Maciejewski, M. Przybylska, Differential Galois theory and Integrability. Internat. 
J. Geom. Methods in Modern Phys., v. 6(8), pp. 1357-1390, 2009. 

[14] I. Marquette, P. Winternitz, Polynomial Poisson algebras for classical superintegrable 
systems with a third-order integral of motion, /. Math. Phys. v.48, p. 012902, 2007. 

[15] J. J. Morales Ruiz. Differential Galois Theory and Non-lntegr ability of Hamiltonian Systems, 
volume 179 of Progress in Mathematics. Birkhauser Verlag, Basel, 1999. 

[16] J. J. Morales-Ruiz, J. P. Ramis, A note on the non-integrability of some Hamiltonian 
systems with a homogeneous potential. Methods Appl. Anal, v. 8(1), pp. 113-120, 2001. 

[17] K. Nakagawa, Direct construction of polynomial first integrals for Hamiltonian systems with 
a two-dimensional homogeneous polynomial potential, Ph.D. thesis. The Graduate Univer- 
sity for Advanced Studies, Japan, 2002. 

[18] E. G. C. Poole, Introduction to the theory of linear differential equations. Dover Publications 
Inc., New York, 1960 

[19] M. Przybylska, Darboux points and integrability of homogenous Hamiltonian sys- 
tems with three and more degrees of freedom. Regul. Chaotic Dyn., v. 14(2), pp. 263- 
311, 2009. 

[20] M. Przybylska, Darboux points and integrability of homogenous Hamiltonian sys- 
tems with three and more degrees of freedom. Nongeneric cases. Regul. Chaotic Dyn., 
V. 14(3), pp. 349-388, 2009. 



36 



[21] M.R Ranada, M. Santander, Complex Euclidean super-integrable potentials, potentials 
of Drach, and potential of Holt, Phys. Lett. A, v. 278, pp. 271-279, 2001. 

[22] A. V. Tsiganov, Towards a classification of natural integrable systems, Regul. Chaotic 
Dyn., V. 11(3), pp. 343-362, 2006. 

[23] A.V. Tsiganov, On maximally superintegrable systems, Regul. Chaotic Dyn., v. 13, pp. 
178-190, 2008. 

[24] A.V. Tsiganov, Addition theorem and the Drach superintegrable systems, /. Phys. A: 
Math. Theor., v. 41(33), p. 335204 (16pp), 2008. 

[25] A.V. Tsiganov, Leonard Euler: addition theorems and superintegrable systems, Regul. 
Chaotic Dyn., v. 14(3), pp. 389-406, 2009. 

[26] A.V. Tsiganov, On bi-integrable natural Hamiltonian systems on the Riemannian man- 
ifolds, .ar)(ivj2Q063914^, accepted to Journal of Nonlinear Mathematical Physics, 2010. 



37 



